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Abstract 

A deterministic computational procedure for describing the transport of electrons in condensed media is formulated to 
simulate the effects and exposures from spectral distributions typical of electrons trapped in planetary magnetic fields. 
The primary purpose for developing the procedure is to provide a means of rapidly performing numerous repetitive 
transport calculations essential for electron radiation exposure assessments for complex space structures. The present 
code utilizes well-established theoretical representations to describe the relevant interactions and transport processes. 
A combined mean free path and average trajectory approach is used in the transport formalism. For typical space 
environment spectra, several favorable comparisons with Monte Carlo calculations are made which have indicated that 
accuracy is not compromised at the expense of the computational speed. 
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1. Introduction 

The concept of developing a rapid analysis electron 
transport code at NASA Langley Research Center (LaRC) 
arose from a desire to have a companion code for the 
LaRC deterministic transport code HZETRN [T|. The 
algorithm developed here for electron transport is based 
on the Continuous Slowing Down Approximation (CSDA) 
combined with an elastic multiple scattering formulation 
to define an electron mean free path and a transmission 
function at a given target location. The cross sections 
relevant to CSDA and multiple scattering are described 
in the following section. The pertinent cross sections are 
calculated with well-established theoretical formulas, valid 
over the energy range representative of trapped electrons 
in planetary magnetic fields. Consequently, a calculation 
of positron production and annihilation is incorporated. 
The cross section equations are cast in terms of specific 
elemental atomic species, from which cross sections ap- 
plicable to any user defined molecular system (mono or 
polyatomic) are automatically constructed. 

Details of the transport formulation are described in 
a subsequent section, and are based on motion related 
to some initial direction (axis of propagation). Behavior 
of slowing and stopping of electrons and their associated 
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bremsstrahlung is evaluated in terms of quantities con- 
tained in the direction of the hemisphere centered on the 
initial direction. The scattering and mean free path cal- 
culations permit definition of mean trajectories relative to 
the initial motion direction, so that axisymmetric spatial 
and energy distributions may be inferred. 

Comparisons with the corresponding Monte Carlo (MC) 
calculations are shown. The MC cases utilized several 
hours of machine time, whereas the LaRC calculations 
were practically instantaneous. Despite some difference 
in final magnitudes, general functional behavior is consis- 
tent in the results. Consequently, the deterministic code 
can be reliably used in trade studies where rapid analyses 
are necessary. 

The present code has the desirable feature that the re- 
quired inputs are user defined through hies specifying ma- 
terial composition along with the incident environment dif- 
ferential energy spectra for electrons and photons. Also 
included with the code is a permanent atomic database 
Hie which details the pertinent atomic properties. Most 
results in this paper are presented for three selected ma- 
terials: aluminum (Al), tantalum (Ta), and water (H 2 O). 
This selection of target materials was chosen to provide 
examples of light and heavy materials in addition to exam- 
ples of elemental and compound materials. The materials 
were also selected for their uses in space applications. Alu- 
minum is a common spacecraft structural material, water 
closely approximates human tissue, and tantalum is used 
in shielding sensitive electronics from electron and photon 
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exposure in the planetary trapped radiation belts. 

2. Cross Sections 

The theoretical description of propagation of ener- 
getic electrons and photons in condensed media requires 
some representation of interaction cross sections with con- 
stituent atoms of the medium. Application of principles of 
quantum and semi-classical physics have gone far toward 
providing formulations for electron and photon interaction 
cross sections. Although many of these formulations are 
mathematically complicated, they are often amenable to 
approximation and parameterizations that greatly simplify 
their practical application. The following subsections de- 
scribe the formulations used here for interactions of elec- 
trons and photons. 

2.1. Electrons 

Free energetic electrons passing through a material are 
slowed and may eventually be stopped as they interact 
with the fields of the nuclei and bound electrons of the 
material. Interactions also result in a change of direc- 
tion (scattering). In the present analysis, the electron 
deceleration process is expressed in terms of energy loss 
by collisions and the accompanying photon production 
(bremsstrahlung). The dominant scattering process is as- 
sumed to be that of small angle elastic (Rutherford) scat- 
tering and is the only scattering process taken into account 
in the present formulation. 

2.1.1. Collisional Slowing and Stopping 

Electron deceleration by collision processes, in which en- 
ergy from the projectile electron is imparted to the elec- 
trons of the medium, is usually cast in terms of the stop- 
ping power, S , which is defined as the energy loss per unit 
scaled distance, — dE/dx . Note that the distance x is mea- 
sured in units of g/cm 2 . The resultant expression for the 
collisional stopping power has been often referred to as 
the modified Bethe-Bloch formula and is given here in the 
notation of Anderson [5j, 

S col = ~ d ^ = 2nr 2 e m e c 2 P \N 0 Z/A)G T C ± , (1) 

where 


and 


{T + m e c 2 ) 2 
T(T + 2 m e c 2 ) 


distance measured in g/cm 2 , N 0 is Avogadro’s constant, Z 
is the atomic number of the material, A is the atomic mass 
of the material, T is the kinetic energy of the electron or 
positron, I is the mean atomic ionization potential, <5 is 
the density correction term, and A is the shell correction 
term. The F ± quantity in Eq. ([2]) is a function only of 
projectile energy and is slightly different (a few percent) 
for electrons and positrons. The present formulation im- 
plements only the F ± function that applies to electrons. 
It may be written as the sum of three terms [2], 
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where 
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The quantity S in Eq. <[2j) represents a modification to 
the Bethe-Bloch stopping power usually referred to as the 
“density correction.” For this work, a parameterized form 
of 5 based on extensive experimental data j3j is used. 

The remaining correction term, A, in Eq. ([2]), is called 
the “shell correction” and is only important at low energies 
when the projectile electron speed, v , approaches the speed 
of a bound electron, vj. Its value, as used in the present 
calculations, is given by 

\~( Vl \ 2 ~ 1(1 + 2m e c 2 )(T + m e c 2 ) 2 

VW ~ T{T + 2m e c 2 )(I + m e c 2 ) 2 ' U 

It is seen in Eq. 0 that S co \ is unbounded as T — > 0 
due to the Gt factor. To correct the unbounded behavior 
as T — > 0, the non-relativistic stopping power derivation 
provided by Bohr and discussed in the paper by Sigmund 
|4j can be used. Bohr’s derivation results in a stopping 
power that approaches zero at very low energies and ex- 
hibits a maximum. The Bohr formula may be written in 
CGS (centimeters-grams-seconds) units as 
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The superscript (±) indicates applicability to both elec- 
trons (-) and positrons (+). Here, r e is the classical elec- 
tron radius, m e is the rest mass of the electron, c is the 
speed of light in vacuum, p' is the scaled mass density with 


with n = pNo I A is the number density, p is the mass den- 
sity, e is the electric charge, C 1 = 2e -7 , 7 is Euler’s con- 
stant, wj = I/h 7 and h is the reduced Planck’s constant. 
The kinetic energy at which the Bohr formula for stopping 
power exhibits a maximum is found from dS ^ hK = 0, 
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From Eq. 
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where Tg is the projectile kinetic energy at peak stopping 
power and exp is the base of natural logarithms. If Eq. 
([I]) is compared with Eq. ([5]), the natural logarithm term 
in Eq. ([H]) can be identified with the C ± in Eq. ([l]). This 
gives 


C t 



and 


-Smax = 2irrlm e c 2 N 0 G T (Z / A)C± ax . 

While Eq. (|T|) applies to T > Tg, an extension of Bohr’s 
stopping power formula to low energies, for T < Tg, is 
provided by Sigmund [4]. For the case of T < Tg, a simple 
functional fit for the collisional stopping power in units of 
MeV cm 2 /g has been derived for the present work, 

Scd = Smax exp {-0.179 [In (T s /T)] 2 05 }. (6) 

Three simple materials important to space exploration 
applications have been selected to demonstrate results 
of the present model calculations: aluminum, tantalum, 
and water. These materials were chosen to represent a 
broad range of atomic properties and material densities. 
The collision stopping powers, as calculated for these sub- 
stances over the energy range of the code formulation, 
were compared with the stopping powers obtained from 
the database tabulations of the National Institute of Stan- 
dards and Technology (NIST) [5]. The comparisons indi- 
cated that the present formulation is quite acceptable over 
the energy range of interest for space applications (0.01 - 
1000 MeV), with the largest discrepancies occurring at the 
lowest energies for heavy elements [6i. 


2.1.2. Radiative Energy Loss (Bremsstrahlung) 

Accelerating (and decelerating) charged particles lose 
energy by photon emission. Elaborate quantum-theoretical 
calculations have gone far toward quantifying these 
bremsstrahlung cross sections. A detailed description of 
the theoretical results is given in Koch and Motz [7] . The 
formulation may be greatly simplified by using several pa- 
rameterizations; the implementation for this work is essen- 
tially the same as that for the default cross sections of the 
EGSnrc Monte Carlo code [8J. Several pertinent param- 
eters that lead to the ultimate cross section formula may 
be defined as [8] 
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Here, E v is the emitted photon energy and a is the fine 
structure constant. The term F c serves to provide a high 
energy Coulomb correction. Another parameter defines 
the effects of screening in terms of an energy variable [8], 
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which involves the electron total energy, E e = T + m e 
The screening functions have been specified as [8 a 
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20.867 - 3. 2425 s + 0.625<5 2 

5 S < 1 

21.12 - 4.184 ln($ s +0.952) 

5 S > 1 

20.209 — 1.935d s + 0.086<5 2 

< 1 
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5 S > 1 


To account for the triplet production process, a param- 
eter, £, is introduced which depends upon atomic charge 
and involves the Coulomb correction above [8j. A param- 
eterized fit for £ in units of cm 2 /g has been developed for 
this work as 

3 

£« 1.147 sin(0.42 In Z) +0.12 |[sin(0.761n Z)] 2 j 4 . (10) 

The final expression for the cross section for production 
of a photon of energy E v by interaction of an electron of 
total energy E e on an atom of charge Z may be written as 


dqy, e _ ar 2 e N 0 Z(Z + £) 
dE u E u A 

where 


(11) 


H v = 


1 + 


E e 


2 ( E e — E v 
E e , 


(0i - F c ) 
(02 - F c ). 


An obvious energy constraint is that only photons can 
be produced that have energy less than the initial electron 
kinetic energy. Consequently, da^^ e /dE v is calculated at 
each electron energy T, for selected values of E„/T: [0.001, 
0.05, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 0.9, 0.95, 0.98, 0.99, 1.0]. 
This distribution of energy values has been chosen to ac- 
curately resolve the cross section variations at both ends 
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of the energy spectrum. A source term £, for production 
of photons by an electron flux ip e , may be written in terms 
of the bremsstrahlung cross section as 


T 

C {E v , T) = J y e {T') d -^dT' for E„ < T. (12) 

E v 

The corresponding energy loss term due to radiative pro- 
cesses may be evaluated as 


*Srad — 


'k^k- 


(13) 


Using Eq. (13 1, the radiative loss stopping powers for 
the three materials selected for this study were calcu- 
lated. In contrast to the collision stopping powers, the 
bremsstrahlung stopping power is greater for higher charge 
elements and increases monotonically with energy. The 
logarithmic stopping power for the present approximate 
formulation is nearly linear with the logarithm of the en- 
ergy and begins to have substantial effects for kinetic en- 
ergies greater than a few MeV. Comparison with NIST 
[5] data revealed some disparities at low kinetic energies 
(T < 0.1 MeV) but, as will be shown subsequently in 
Sec. |2.3[ bremsstrahlung effects on the general transport 
process have a significant impact only for kinetic energies 
greater than several MeV. 


where 0 is the scattering angle and e is a parameterized 
screening function, 

e « 1413[(1 + T/m e c 2 ) 2 - 1 }{Za)~ 2/3 . 

Integrating Eq. ( fl5] ) over all scattering angles to obtain 
the total cross section allows the electron transport mean 
free path to be written as |Bj 


Atr — p / no - r 

A (T 2 + 2Tm e c 2 ) 2 

2t rIV 0 e 4 Z 2 [ln(2e) - 1] (T + m e c 2 ) 2 ' 1 J 

The above quantity will be dealt with in the subsequent 
section on the transport algorithm. Note that the behav- 
ior of the transport mean free path increases monotoni- 
cally with kinetic energy and is nearly logarithmically lin- 
ear with the log of the kinetic energy, having generally 
higher values for lower atomic weight elements. 

Another parameter related to elastic multiple scattering 
is the mass scattering power m- This quantity is analo- 
gous to stopping power, but refers to solid angle scattering 
rather than energy loss. The mass scattering power is de- 
fined in Li and Rodgers m as 

Nn 

-fy{0r), (17) 

where 


P, = TV 


2r e m 2 c l Z(T + m e c ) 
T(T + 2m e c 2 ) 


2.1.3. Multiple elastic scattering 

Strict conservation principles require that electron-atom 
interactions generally result in energy exchange accom- 
panied by re-direction with respect to spatial variables. 
Approximations that tend to decouple energy loss pro- 
cesses and directional changes greatly simplify electron 
transport analysis. In the present work, these processes 
are made practically independent. Energy loss is specified 
by collision and radiative losses without regard for direc- 
tional change, while projectile trajectories are described by 
elastic scattering interactions. Such assumptions call for 
careful scrutiny with regard to both broad energy spectra 
and the wide variety of material types. Some precedents 
have been set in earlier works [91 [10] in which such ap- 
proaches have been used. In the present formulation, the 
elastic (Rutherford) scattering cross sections, or, are im- 
plemented as a basis for calculating an electron transport 
mean free path, 

Atr = (14) 

na R 

where n = pN^/A is the number of scattering centers per 
unit volume in the medium. The present development 
closely follows that of Kawrakow and Rogers [5] , in which 
the screened Rutherford cross section is given as 

do-R _ 0 T Ze 2 (T + m e c 2 ) l 2 

d9 ~ * [T(T + 2m e c 2 )(l-cosd + l/e)\ ’ 1 ’ 


Y{6 r ) = ln(l + ^)-l +rT# 

with 

_ y/T(T + 2m e c 2 ) 
r am e c 2 Z 1 / 3 

The mass scattering power is used in the transport al- 
gorithm to relate electron traversal on a specific trajectory 
to the axial penetration distance. The above formulation 
predicts a nearly linear monotonically decreasing variation 
of the logarithm of scattering power with the logarithm of 
energy. 

2.2. Photons 

Bremsstrahlung photons produced by electron-atom in- 
teractions are attenuated and absorbed in the transport 
medium. In the present work, only the three most im- 
portant processes involved in photon energy degradation 
and absorption are considered: photoelectric absorption, 
incoherent scattering, and electron-positron pair produc- 
tion. These processes are dominant for different energy 
ranges, but when taken collectively, they closely approx- 
imate the total attenuation for the energy range consid- 
ered here (0.001 - 1000 MeV). The contribution due to 
coherent (or Rayleigh) scattering has not been included in 
this analysis since it does not significantly add to the total 
cross section over the energy range of interest (0.001 - 1000 
MeV). 
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2.2.1. Photoelectric Absorption 

Photoelectric absorption involves interactions of a pho- 
ton with a bound electron in which the photon energy is 
absorbed with a corresponding kinetic energy increase for 
the electron. The detailed evaluation of the photoelectric 
cross sections requires the use of complicated atomic mod- 
els in sophisticated solutions of the Schrodinger equation 
for all energies and elements under consideration. For the 
present formulation, liberal use is made of parametric fit- 
ting formulas devised from the extensive calculations and 
tabulations provided by NIST |5j. First, the energy values 
for the absorption discontinuities of the K and L shells 
have been approximated as 

E k « 5.435 x icr 3 Z 2 - 2038 keV (1 < Z < 92), (18) 

E l « 1.5754 x icr 4 Z 2 - 5984 keV (30 < Z < 92). 

Note that energies in this section are in keV, and A-edge 
energies are only considered for Z > 30. In addition, only 
the discontinuities due to the K and L shells have been 
considered. The discontinuities (edges) due to the elec- 
tronic shells occur at values of the photon energy corre- 
sponding to the binding energies of electrons involved in 
the interaction. When the incoming photon has an energy 
greater than the binding energy of a given atomic shell or 
subshell, a discontinuity may occur. 

At the edge discontinuities, the NIST data has also been 
used to develop formulas for the upper (w) and lower ( l ) 
values of the photoelectric cross section: 


a Ku 

= 7165.5£^ 1 - 5474 

cm 2 / g, 

(19) 

a ki 

= 593.59A^ 1327 

cm 2 / g, 


a Lu 

= 9620.6A)" 1 - 5113 

cm 2 / g, 


au 

= 1633.5A)" 1 - 2124 

cm 2 /g. 



Examination of typical photoelectric cross sections re- 
veals that the logarithmic slope, 

_ d(ln ope) 
m “ d(ln E v ) ’ 

is very nearly constant for a given atomic system between 
absorption edges. Using NIST [5] cross section data, 
parametric fits were used to approximate the logarithmic 
slopes. 


For E v > E K - 


m K 


6.9068 x 10~ 3 Z — 3.014 Z < 10 
-3.3147Z- 0 08156 Z > 11 


For .1 keV < E v < Ek' 


m L = 



(21) 


For all Z and E k < E v < 10000 keV: 


<7 PE = Oku 



For Z >29 and El < E v < Ek'- 


ctpe = CTlu 



( 22 ) 


(23) 


An extension to energies less than the lowest edge en- 
ergy considered was extrapolated with a power law func- 
tion approaching an asymptote based on an extrapolated 
maximum absorption cross section. For Z < 29 and for 
energies below the A'-edge, an asymptotic cross section is 
calculated using the average ionization potential, 


f I \ mK 

a] - aKl VFhJ ' 

The extrapolation formula is then 

< ’™ = 1 + A P ES’ 

where the exponent p and coefficient Ap are given by 


(24) 


(25) 


p= (™M)-i’ Ap= < ' <26) 

These constants are found by matching the cross section 
and the slope at the AT-edge. For Z > 29, the low energy 
extrapolation below the A-edge is found by using an L 


subscript instead of K in Eqs. ( 22 )-( 26 ) . 


For the three selected materials, the photoelectric cross 
sections, as calculated with the present parameterization 
were compared with their NIST [5] counterparts and agree- 
ments were considered to be fair to good [5]. 


2.2.2. Inelastic Photon Interactions 
As the photoelectric absorption coefficient decreases 
rapidly with increasing energy, the next process to domi- 
nate is when the photon imparts a portion of its energy to 
an electron resulting in scattering of a lower energy photon. 
The process is referred to as inelastic (incoherent) scatter- 
ing. The cross section formula for release of an electron of 
kinetic energy T e by a photon of initial energy E u has been 
taken from Anderson [2] as the free electron Klein-Nishina 
process. The energy constraints are 0 <T e < T max where 

rr 2 El 

max “ m e c 2 + 2E V ' 

Here, T e is the kinetic energy of the emitted electron and 
A„ is the incident photon energy. The cross section for- 
mulation is 
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In order to render the cross section applicable to bound 
atomic electrons, elaborate quantum theory calculations 
are required. Such considerations are only appreciable at 
lower photon energies. For this work, a parametric mul- 
tiplying function has been developed to approximate this 
effect, 


In $(Z,E V ) 


Bi In (Z) + B 2 
E v 


1 [B3-B4 In (Z)] 


(28) 


where B t = 0.471, B 2 = 1.5184, B 3 = 1.003, and 
B 4 = 0.13317 when E v is given in units of MeV. The final 
inelastic scattering cross section then becomes: 


da m 
dT e 


$>{Z,E V ) 


d(JKN 
dT e ' 


(29) 


For each E v , the total inelastic cross section is found by 
integration over the appropriate emitted electron energy; 


Qi = [E 2 + + (E„ - E+) 2 } [0i- F c ], 


and 


Q 2 = 2 -E + {E v -E + )[<j> 2 -F c }. 

The terms £, </>i (S'), ^ 2 ( 6 '), and F c are the same as for 
the bremsstrahlung cross section, Eqs. ([ 7 |)-([To|, with the 
exception that the value of S' is determined as 


136m e c 2 E v 
Z 1 / 3 E + {E V — E + ) ' 


(32) 


Energy conservation constrains the positron energy as 


TH e c t if 1. y r max T / 1 1 p c , 


with 


T — E — c 2 

± max . 


En;ix 

°"'= / d W dr - (30) 
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As a function of incident photon energy, the total inelas- 
tic cross sections, calculated with the present model, were 
compared with the corresponding NIST [F values. The 
model approximation was in better agreement for lighter 
materials. An exception occurred at higher energies where 
the model formulas are slightly inaccurate [8]. However, 
in this energy range, the inelastic cross section does not 
contribute appreciably to the total photon attenuation, as 
will be subsequently shown in Sec. |2.3| 


2.2.3. Pair Production 

Electron-positron pair production may occur when a 
photon of sufficient energy interacts with a strong localized 
Coulomb field. The photon initial energy is transformed 
into the combined rest mass of the newly created particles 
(2 m e c 2 ), with any remaining energy appearing as kinetic 
energy of the new particles. A comprehensive description 
of the process may be found in Motz et al. [T2] , where a 
variety of cross section quantities are derived for several 
aspects of the phenomenon. The formulation from Li and 
Rogers El. chosen for use in the present work, applies to 
the cross section for production of a positron of total en- 
ergy E + by a photon of initial energy E v . Only positron 
production associated with pair production is considered 
in this work. The positron production differential cross 
section may be written in simple terms as 


where 


d(T pp 

~dE + 


N 0 ar 2 e Z(Z + Q 
A El 


0 Qi + Q 2), 


(31) 


The above formulas for pair production become some- 
what inaccurate as the pair production threshold is ap- 
proached. Consequently, a power law function is fit at the 
first energy grid value greater than 4 MeV (= E 4 , see be- 
low). Now, the low energy pair production cross section is 
given as 


(r pp = A pp (E y - 1.03) y for 1.03 MeV < E v < E 4 , (33) 


with 


dcr. 


V = 


pp 


dE, 


E+=Ea 


E a - 1.03 

0-pp(^4) 


and 


(£4) 


A = W PP 

pp (E 4 -1.03)v' 


A source term for positron production may be derived 
in a manner similar to that described for production 
of bremsstrahlung photons. See Eq. (12). For a pho- 


ton differential flux spectrum ip v (E v ,x), positron produc- 


tion may occur between energies E p 


= 2m e r and 


A'pp^max — Ei/, max 


— 2m e c , where E v , max is the highest 
photon energy value considered. The expression for the 
source term becomes 


^pp,max 

C(E+,x)= j ip v {E' v ,x)^-dE' v . (34) 

E 

pp,min 

It should be noted that in parallel with the positron pro- 
duction process, a similar cross section expression applies 
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to the creation of the partner electron, along with the cor- 
responding source term for secondary electrons. Just as 
for the inelastic process, the pair production cross section 


may be found by integrating Eq. (31 1 over the allowed 
positron energy ranges, namely from m e c 2 to E v — m e c 2 . 

Results of pair production cross sections, as calculated 
from the present model were compared with the appro- 
priate NIST [S] values. Agreement was very good except 
for energies near the threshold region where the power law 
extrapolation formula was invoked [B]. 


2.3. Cross Section Implementation 

The individual cross sections described previously are 
used in combination to provide three critical parameters 
essential for description of the electron-photon transport 
process. The ultimate slowing and stopping of electrons is 
governed by the total stopping power obtained from the 
sum of the collision and radiative stopping powers. The 
collision process dominates at low energies (less than ap- 
proximately 1 MeV). Here, the radiative process assumes 
the dominant role at high energies (greater than approx- 
imately 1 MeV). A transitional minimum occurs in the 
neighborhood of 1 MeV where the lower energy inaccura- 
cies in the radiative stopping power of the present model 
have insignificant impact on the total stopping power. 

The general attenuation and extinction of photons in 
a medium is found to be closely related to the sum of 
the photoelectric, inelastic, and pair production cross sec- 
tions. In the present model, the coherent (Rayleigh) scat- 
tering processes were not considered because of their rela- 
tive unimportance in space radiation effects. The inaccu- 
racies in the inelastic cross sections at energies larger than 
approximately 100 MeV have insignificant impact on the 
total attenuation of photons due to the pair production 
process dominating at higher energies jB]. 

The final parameter of importance to the electron- 
photon transport process is that of the photon energy de- 
position coefficient which is used to calculate effective pho- 


ton dose. See Eq. (63) in Sec. 4.1 This relates to those 


fractions of the photon attenuation processes (/pe> /in, 
/ pp ) that produce secondary electrons. For photoelectric 
absorption, it is assumed that the secondary electron is 
emitted with the same energy as the incident photon (i.e. 
the binding energy of the released electron is neglected) 
and the value of /pe is unity. In the case of inelastic scat- 
tering m , 


where E„i n is the energy of the photon produced in the 
inelastic process initiated by a photon of energy E„. For 
the pair production process, / pp = 1 — 2 m e c 2 /E„. 

An additional reduction in the photon energy deposition 
coefficient is required to discount secondary electrons that 
subsequently produce additional bremsstrahlung. In the 
present model, this factor is represented by (1 — SVad/STot), 


where St ot = S ra( j + S co i . This factor is only applied to 
the inelastic and pair production processes, so that the 
model expression for the total photon energy deposition 
coefficient becomes m 

T 

max, in 

Me„=-PE+ J 

0 

Tmax,pp 

+ / a-^/pp^r, (35) 

o 

where the upper limits on the integrations apply to the 
maximum allowed emitted electron energies for the re- 
spective processes. The energy deposition coefficient is 
described in detail in the X-Ray Data section of the NIST 
website [5]. 


3. Transport Formulation 

In the present work, the general transport process is 
focused principally on the description of the penetra- 
tion of a primary electron field along with secondary 
bremsstrahlung photons generated by electron-atom 
interactions. For electrons, an essentially one-dimensional 
formulation is developed with some reference given to 
the effects of radially symmetric scattering. While the 
electron propagation is determined by the stopping 
powers, the photons are assumed to be transported along 
the direction of travel of the electrons. Photon intensity 
is governed by the radiative transfer equation El, which 
utilizes the calculated attenuation coefficients. 


3.1. Electrons 

A quantity of fundamental importance in describing the 
transport of electrons in matter is the maximum distance 
of travel, which is determined from the energy loss stop- 
ping powers calculated from the modified Bethe-Bloch and 
bremsstrahlung formulas in Eqs. 0 and ( |13[ ). This quan- 
tity is referred to as the CSDA range R(E). For a given 
energy, E , the CSDA range is defined as 


R(Ji)s /£ = /ra? (36) 

o o 

The inverse of the R(E) function, R~ l (E) = E(R), may 
be used to construct the variation of dE/dx as a function 
of distance traveled by an electron of initial energy Eq. 
Such a function, S(Eq , x), is analogous to the usual Bragg 
curve, but neglects straggling and measures energy loss as 
a function of distance traveled along the particle’s path 
rather than perpendicular depth in a material. Straggling 
is accounted for after the incorporation of the statistics 
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associated with scattering. The residual energy of an elec- 
tron, at position s, is given by 


R 


W(s) = j S(Eo, x)dx — j S(Eo,x)dx 
o o 

S 

= E 0 - J S(E 0 , x)dx, s < R. (37) 


For a beam of mono-energetic electrons, it is assumed 
that the maximum distance traversed is the CSDA range 
and straggling is neglected. However, the average distance 
of penetration along the beam axis direction is generally 
less than the CSDA range due to multiple scattering ef- 
fects. The average deflection for a unit path length can 
be defined by the inverse mean free path (A)) 1 ) [5] which 
yields 


and 



. o 


(38) 


S 



0 





0 


(39) 


In general, (x) is defined to be the arithmetic average of 
the variable x. For electrons, ( z(R )) is interpreted as an 
“effective range” or the axial distance at which 50% of the 
electrons have stopped and takes place at the z-value for 
which the transmission is 0.5. 

A further critical assumption is made by specifying a 
Gaussian distribution about ( z{R )), where ( z(R )) is the 
distance of peak electron population. In addition, it is 
assumed that practically all electrons are stopped at z = 
R(E). This Gaussian represents the variation in electron 
path length due to multiple scattering effects. By invoking 
the formula for half-width of the Gaussian distribution, a 
value for deviation, H, is found as m 


H , R-{z{R)) 

2^2 In 2 


(40) 


The resultant Gaussian function may be interpreted as the 
probability that an electron penetrates an axial distance 2 
within increment dz, 


p , v i / 

P{z ) = ex P i - 




2tt 2 


(41) 


Integration of the Gaussian function results in an er- 
ror function prescribing the number of particles having 


stopped over distance z. The complementary error func- 
tion [16j results in the corresponding transmission func- 
tion, 


rj(z) = 0.5 < 1 — erf 


z - (z(R)) 

H^2 


(42) 


For simplicity, a special algorithm for the standard error 
function has been utilized in the present code Ba- 
in addition to specifying the transmission probability of 
electrons as a function of axial penetration distance z, it 
is necessary to evaluate the variation of primary electron 
energy as a function of z. A formula has been developed 
m in terms of the moments of the mass scattering power 
that expresses, in cylindrical coordinates (z,r), the dis- 
tance covered for an average electron trajectory having 
axial distance z and deflection radius r. The moments are 
calculated as 


z 

= J P s {u)(z - u) z du, 


(43) 


with P s given by Eq. (17 1. The moments are then used 
to derive an equation |T(7 for distance along an average 


trajectory passing through (z,r), 


Z 

s(z,r) = j 


1 + ao(z ) + 


/ r 2 
\a 2 (z) 


- 1 


qj(z') 

02(2') 


dz . (44) 


The quantity yj a 2 {z) has been identified [TD] as the half- 
width at half-height of the distribution of deflection radius 
at axial distance z. This observation leads to the infer- 
ence that a “most typical” trajectory intersects the point 


(Zjy/cRz). The third term in Eq. (44) then vanishes. An 
“average mean path” now becomes, 


z 

{ s{z )) = J \Jl +a 0 {z')dz 


(45) 


Reference to Eqs. (39) and (45) shows that the relation- 


ship between the mean axial distance and average mean 
path is such that (z) < z < ( s ). An average residual 
energy {W) may be associated with an average distance 
along a path length ( s ) by again interpolating on a table of 
CSDA range-energy relations. However, the interpolated 
value ( W ) reaches zero before z = R is attained because 
( s ) is a value that includes average deflection and z = R is 
the axial distance traveled in the continuous slowing down 
approximation assuming no angular deflection. The axial 
penetration, however, may extend beyond the distance at 
which ( W ) vanishes due to energy straggling. A provision 
for energy straggling is supplied by introducing a power 
law function of the form 


(W)=A w (R-z) k 


(46) 
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This formula is invoked somewhat arbitrarily for 2 -values 
for which (W) < 0.5 VFcsda- When this condition per- 
tains, the constants A w and k are evaluated by point-slope 
matching at the 2 -value where (W) = 0.5 VFcsda- 

In going from the boundary to distance z , particle num- 
ber conservation requires that ip(0,T o )ST o = ip(z,T)ST, 
where the energy scaling corresponds to that prescribed 
by the CSDA range-energy relation for the material. In 
general, SY is defined as a small, finite difference in the 
value of Y. The above transport parameters, along with 
particle number conservation, may be used to devise an 
expression for the differential flux spectrum of electrons, 
<p(z,T), at distance 2 with kinetic energy T, given the 
initial spectrum at z = 0 with initial kinetic energy To, 
t( 0,T o ): 


tp(z,T) = 


tp(P,T 0 )r,(z)S(To) 

S(T) 


(47) 


In addition, the transmission function r/(z) is used to ac- 
count for particles not arriving at 2 due to interaction not 
accounted for in the stopping power. Division of the ki- 
netic energy increments by distance increments leads to an 


expression in terms of stopping powers in Eq. ( 47 ) . Thus 


for a slab of thickness w, the differential spectrum of pri- 
mary electrons may be found at any axial distance z < w. 

If primary electrons are the only concern, the calcula- 
tion could end with Eq. (47). However, bremsstrahlung 


photons are generated during the transport process. In 
order to take this process into account, a spatial grid is 
established within the slab layer. In the present code, a 
slab layer is assigned a spatial grid, normally of 20 to 30 
points, with spacing increasing monotonically according to 
the following formula, 


dt = 


i — 1 

N- 1 


(48) 


Here, di is the depth in the material slab of the ith point 
and N is the total number of points in the spatial grid of 
the slab. The photon source term may then be calculated 
at each spatial grid point as 

T(di) 

C{E v ,d i )=J<p e (T')^-dJ / . (49) 

E v 

Finally, the primary electron spectra provided by the 
present code are calculated at each spatial grid point as 


VeMiT) = 


y(0, E 0 )r](di, E 0 )S(E 0 ) 
S(T) 


(50) 


Implementation of Eq. ( 47 ) for the illustrative low Earth 


orbit (LEO) spectrum [181 normally incident on aluminum 
provides the results shown in Fig. [l] (the environment spec- 
trum is labeled Z = 0.0). The stopping of low energy elec- 
trons and the transfer of high energy electrons to lower 
energies is evident in the spectral functions at increas- 
ing penetration depths. Similar results for Jovian trapped 



Figure 1: LEO electron spectra at several aluminum thicknesses. 


electrons in the vicinity of Europa in aluminum are shown 
in Fig. fusing an environment from the NASA Jet Propul- 
sion Laboratory (JPL) Galileo Interim Radiation Electron 
(GIRE) Model [T5]. The much greater penetration and 
more gradual attenuation is clearly evident for this very 
high energy environment. The peak in the flux occurs be- 
cause low energy electrons are more likely to be stopped in 
aluminum which appears as a decrease in flux at low en- 
ergies. High energy electrons, on the other hand, will be 
shifted to lower energies but not lost from the spectrum 
which leads the shift in the peak of the flux to higher en- 
ergies with a decrease in peak flux. 

3.2. Photons 

In the previous section on cross section evaluation, it 
was shown that a spectral source term for bremsstrahlung 
photons may be calculated. See Eq. ©• When the dif- 
ferential electron flux spectra is specified as a function of 



Figure 2: Europa electron spectra at several aluminum thick- 
nesses. 
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axial distance z, a corresponding photon flux may be de- 
termined using a solution of the radiative transfer equa- 
tion for an emitting-absorbing medium M- In addition 
to the source function, the total attenuation coefficient is 
required. The total attenuation coefficient is the sum of 
the relavant cross sections, 

OTot = °PE + Om + CTpp. (51) 

The photon flux spectrum at each spatial grid point di 
then involves integration of the source term multiplied by 
the attenuation factor; 


and t representing the positron kinetic energy in rest mass 
units (0.511 MeV). 

Positrons initially at spatial location Xi arriving at Xi + 
S x with energy Ej began traversal of this increment at a 
higher energy, E 0 = E(R + Sx), where R is the residual 
range at a \ + fa. The CSDA process over the interval Sx, 
with no loss of positrons requires that 

tp(E 0 ,Xi)S(E 0 ) = (p(Ej,Xi + i)S(Ej). (54) 

When annihilation takes place, the fraction of positrons 
from the incident spectrum lost in distance Sx is 


di 

Vv,i(d u E v ) = J £(x, E v ) exp [-a To t(di - x)] dx. (52) 
o 

Photon spectra for the case of the Europa spectrum are 
illustrated in Fig. [3] and correspond to the electron spec- 
tra Fig. [2] The persistence of higher energy photons at 
large aluminum depths is indicative of their highly pen- 
etrating nature. The minimum and maximum values of 
the flux at low energies are due to the discontinuity in the 
photoelectric absorption cross section. 


3.3. Positrons 


Ordinarily, exposure effects due to space environment 
electrons may be adequately evaluated from specified elec- 
tron and photon flux values at a given target location. 
The pair production contribution to the energy absorp- 
tion coefficient should provide a fair estimate of exposure 
due to positrons. However, if details of specific positron 
interactions are of interest, it becomes necessary to include 
explicit transport of the positrons and their associated an- 
nihilation photons. In the present code, source terms and 
flux for positrons are evaluated, along with source terms 
for the annihilation photons, in an uncoupled manner with 
a view toward allowing the user to assess the relative im- 
portance of the presence of positrons. The positron source 


term, £(E + ,x), has already been given as Eq. (34). These 


antimatter particles have a “sink” term by annihilation 
with the constituent bound electrons of the medium. This 
term may be represented as [5] 


A(s cr ann (r) 


N e 7T r 2 
t + 2 


t 2 + 6r + 6 
t(t + 2) 



(53) 


where N e is the number of electrons per unit mass of ma- 
terial and (j a nn ( t ) is the total annihilation cross section, 
with 


h\ = In 


T + 1 + yj t(t + 2) 


hi 


T + 4 


r y 


1 

1 / Al e <7 aIm 

E ^Rj + x ) 

dx 

l 0 


J 


The positrons transmitted through interval Sx is then 


( Ej,Xi + Sx) = psAt{8x)p [E ( Rj + Sx) , a;*] , (56) 

with 

S [E ( Rj + fa)] 

PS= S(E~) ' 

Additional contributions to the positron flux at x , + 1 = 
Xi + Sx arise from the distributed sources in the interval 
that must include an additional annihilation factor, As; 


( Ej,x + Sx) 


with 


Sx 

J £[E(R+ x ),Xi+x ]As(a: )dx , (57) 
o 


A s (x ) = exp { - / N e a a 


E ^ Rj + x ) 


dx 


(58) 
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Figure 3: Photon spectra from Europa electrons m at several 
aluminum thicknesses. 



Kinetic Energy, MeV 



Energy, MeV 


Figure 4: Positron flux at selected thicknesses in aluminum for 
Europa spectrum m ■ The curves for a depth of 8.89 g/cm 2 
and 15.3 g/cm 2 overlap above 1 MeV. 


Figure 5: Annihilation photon source for Europa spectrum m 
in aluminum. 


Here, As is the fraction of positrons lost from the positrons 
created in the distance increment. 

The total positron flux, ip e + = + ip( s \ is calculated 

in the code by numerical integration procedures applied to 
each of the established spatial grid intervals. Results for 
the Europa spectrum for selected thicknesses in aluminum 
are shown in Fig. [4] and may be compared with the corre- 
sponding photon flux of Fig. [3j The peak in the positron 
flux shown in Fig. [4] occurs at a positron kinetic energy 
of approximately 1 MeV. This peak in the positron flux 
is due to a minimum in the stopping power model at 1 
MeV [6]. It is seen that the positron flux is a relatively 
small fraction of the governing photon field and that the 
uncoupled calculation is justified. 

When the positrons are annihilated, photons are gener- 
ated in accordance with the differential cross section [5]. 


Finally, a source term for the annihilation photons per unit 
mass of material may be expressed as 


J- max 

C/,ann {E v ,x) = J ip(E+, x)^^~dT + . 


(61) 


Production of annihilation photons occurs with the 
highest probability near the electron/positron rest mass 
energy. For the Europa case of Fig. [4j the annihilation 
photon source terms calculated according to Eq. (61) are 


shown in Fig. [5] and exhibit the expected peak at 0.511 
MeV. Since this contribution is small relative to the other 
processes, these photons are not transported in the present 
formulation. 


4. Dosimetric Calculations and Comparisons 


dcr ann 
dE v 


nr 2 e 

t(t + 2) 


[q(K) + q(r + 2 


«)] . 


(59) 


where 


s I 


+ 2 + 2 


T + 1 
r + 2 



- 1 , 


with r and n being the positron kinetic energy and the 
photon energy, respectively, in units of electron rest mass. 
The annihilation photons are restricted by kinematics to 
an energy range of 


m e 
l + a' 


< E 


i/,ann 


< 


m e 

1 - a'’ 


(60) 


where 


a 


/ 


r + 2 


The conventional radiation dose for ionizing radiations 
in matter is defined as the energy imparted to a mass 
element of the material by the particles traversing the 
elemental mass. In general, only energetic charged 
particles are responsible for the energy transfer. The 
“dose” attributed to uncharged particles (e.g. photons, 
neutrons) results from the charged particles generated 
within the medium by the neutrals. 


4-1. Electron and Photon Dose Evaluation 

Energy deposition of electrons is calculated by multi- 
plying the local flux (differential in energy) by the total 
stopping power and integrating over energy. The dose at 
a given location is given by 

OO 

D e (x) = J <p e (x,E)S(E)dE. (62) 

0 
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In Eq. (62), E represents the electron energy at position x. 
The local differential flux, ip e , has units of electrons/ (cm 2 
MeV). The stopping power, S, has units of MeV cm 2 /g, 
while the units of D e are MeV / g which may be converted 
to cGy or rads upon multiplication by the conversion fac- 
tor 1.602 x 10 -8 . For photons, absorbed dose is charac- 
terized by the energy deposition coefficient, /x en , which is 
described in Eq. ( 35 1 and is expressed herein as a mass ab- 
sorption coefficient with units cm 2 /g. The dose attributed 
to photons is then 


OO 

^ fjjen Ejy ^p v lEvt x^dEi/. 


(63) 


The use of /r en to calculate effective photon dose (as is 
done in the present code) is a simplistic approximation 
method that has been widely used and considered to be 
an adequate representation of this dose contribution 00. 
This approximation assumes that all the energy deposited 
by photons is deposited by charged particles at the point 
of their production. A more direct evaluation would be 
to use the appropriate cross sections for generation of 
secondary electrons (ctpe, cr in , and er pp ) to compute re- 
spective source terms in the manner expressed in Eq. ( 34 ) . 


The corresponding flux terms could then be determined 
using procedures analogous to that described for the 


positron transport in Eqs. ( 53 )-( 57 ) . Such modifications 


may be considered for future upgrades if the additional 
complexity appears to be warranted. 


4-2. Sample Calculations and Comparisons 

Calculations performed with the present code were se- 
lected to examine the behavior of trapped electron spectra 
interacting with light and heavy materials. The chosen 


electron environment spectra (boundary conditions) are 
shown in Figs. |T] and [2] and represent, respectively, the 
environment in LEO (400 km, 51° inclination) |T8] and 
the environment in the vicinity of the orbit of the Jovian 
moon Europa (19j . The initial environment is shown as the 
curves labeled Z=0.0 in Figs. [T] and [2] along with the flux 
spectra at various depths in aluminum. After transport 
through the shield medium, the doses have been evaluated 
in silicon to simulate exposure of a solid state device. 

The Integrated TIGER Series (ITS) Monte Carlo code 
[Ml was used to calculate dose for the LEO spectrum elec- 
trons on aluminum at normal incidence for several depths. 
The identical scenario was used with the LaRC determin- 
istic code and the results are compared in Fig. [6j The 
LaRC code gives generally lower values for the dose ver- 
sus depth functions than the Monte Carlo results, but the 
functional behavior is very similar. Greater differences are 
to be expected for the low energy spectra of LEO because 
scattering processes are more prominent in the transport 
process and are treated differently in the two calculations. 

For higher energy electron spectra, sample application 
comparisons have been made for a Jovian electron environ- 
ment generated by the NASA-JPL GIRE model (TS|. The 
relevant spectrum is shown in Fig. [2] as the curve labeled 
Z = 0.0 and represents electron flux in Jupiter’s equatorial 
plane at the average orbital distance of Europa. The nor- 
mal incidence on a semi-infinite slab scenario was used for 
the materials aluminum, tantalum, and copper-tungsten 
(50%-50%) alloy. The previous LEO spectrum exhibited 
very few electrons above 5 MeV, whereas the Europa elec- 
trons may have substantial population up to and beyond 
100 MeV. 

The dose versus depth curves for the Europa spectrum 
are given in Figs. 13 i and [9] for the specified materi- 
als. In general, the LaRC deterministic calculations show 
improved agreement with the corresponding Monte Carlo 




Figure 6: Comparison of LaRC results in aluminum with ITS 
|20| for LEO spectrum [18|. The results for ITS and LaRC over- 
lap for electrons above approximately 2 g/cm 2 of aluminum. 
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Figure 7: Comparison of LaRC results in aluminum with ITS 
120] for Europa spectrum m 



Tantalum, g/cm 2 

Figure 8: Comparison of LaRC results in tantalum with ITS 
1201 for Europa spectrum m- 

results compared to the lower energy LEO environment. 
The range of scaled thickness for the Europa cases is from 
0.5 to 20 g/cm 2 . It is seen that electron stopping is prac- 
tically complete at the end of this thickness range, after 
which the bremsstrahlung contribution dominates. 

5. Concluding Remarks 

The electron/photon transport code described here has 
been developed with a view toward simplicity and speed 
in analysis of exposure from space environment electrons. 
In its present form, the code may be implemented to great 
advantage in shield material trade studies, numerical and 
statistical experiments, uncertainty analyses, etc., despite 
some observed differences with Monte Carlo comparisons. 
Numerous comparative calculations, in addition to those 
reported here, have been performed with similar degrees 
of agreement as those presented in the text. In particular, 
calculations for the trapped Jovian environment compare 
more favorably with corresponding Monte Carlo results 
than do the comparisons for the much lower energy-range 
spectra of LEO. This fact suggests that the deterministic 
formulas for low-energy scattering should be improved in 
future versions. 

There are several areas in which improvements to the 
code are being considered. The various cross section com- 
parisons with NIST [5] calculations show several instances 
where the LaRC code cross sections may be improved 
by introducing selective correction terms, even though 
the cross section deviations noted have little impact on 
final exposure results. With regard to transport, the 
present code is amenable to immediate extension to a two- 
dimensional axi-symmetric representation. In addition, 
more explicit details pertaining to very high energy pro- 
cesses involving positrons and their annihilation photons 
may be obtained by a formal coupling of these species to 



W-Cu alloy, g/cm 2 


Figure 9: Comparison of LaRC results in W-Cu alloy with ITS 
}20] for Europa spectrum 1191 . 

the general electron/photon transport. Although no sig- 
nificant errors are indicated as a result of the decoupling of 
positron and annihilation photon transport for the Jovian 
spectral environments, use of the present code for high en- 
ergy beam simulation or cosmic ray shower analysis would 
most likely require full coupling of these processes. It is 
therefore natural to anticipate code upgrades and exten- 
sions in the near future. 
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